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1 INTRODUCTION 



ABSTRACT 

We present a clustering analysis of the Texas-Oxford NVSS Structure (TONS) radio 
galaxy redshift survey. This complete flux-limited survey consists of 268 radio galaxies 
with spectroscopic redshifts in three separate regions of the sky covering a total of 
165 deg^. By going to faint radio flux densities (si.4 ^3 mJy) but imposing relatively 
bright optical limits {E k R Ril9.5), the TONS sample is optimised for looking at 
the clustering properties of low luminosity radio galaxies in a region of moderate (0 
^ z ^ 0.5) redshifts. We use the two point correlation function to determine the 
clustering strength of the combined TONS08 and TONS 12 sub-samples and find a 
clustering strength of ro(2:)=8.7±1.6 Mpc (/i=0.7). If we assume growth of structure 
by linear theory and that the median redshift is 0.3, this corresponds to ro(0)=11.0±2.0 
Mpc which is consistent with the clustering strength of the underlying host galaxies 
2.5L-k ellipticals) of the TONS radio galaxy population. 
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Radio galaxies are ideal probes of large-scale structure as 
they are biased tracers of the underlying mass and can 
be easily detected out to high redshifts. By using biased 
galaxies populations, one can efficiently trace huge super- 
structures (i.e clusters of clusters of galaxies) which are still 
in the linear regime and can therefore be directly traced back 
to rare fluctuations in the initial density field at recombina- 
tion. However, in order to be useful probes, it is vital to 
understand how different populations of radio galaxies trace 
the underlying dark matter (i.e. their bias) and how this has 
evolved with time. 

Different populations of galaxies demonstrate a range 
of clustering strengths as measured by the clustering length, 
ro. For example r p ~ 8.1±0.9 Mpc (/i=0.7) for L* elliptical 
galaxies (Norbcr g et alJl2002^ and r p « 5.4 Mpc for IRAS 
selected galaxies jFisher et al]ll994) . The power-law fit to 
the correlation function can also be applied to galaxy clus- 
ter surveys. For clusters of galaxies selected from the APM 
survey, the correlation length has been measured to be rp 
20 4±3. 4 Mpc for Ab ell Richness Class (ARC) >1 clus- 
ters JPalton et al.lll994l) and rp ^ 30.4 Mpc for ARC >2 
JCroft et al.lll997f) . Hence as o ne would expect from high- 
peak bias effects ^lKaiser^ll984^ ■ clusters are more strongly 
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clustered than galaxies. The clustering of radio galaxies will 
follow that of the underlying host galaxy population but it 
may also be affected by the radio triggering mechanism in 
different populations and/or environments. 

Because of the difficulties in obtaining large, complete 
radio galaxy redshift surveys, the clustering properties of 
radio galaxies are relatively poorly constr ained. For local 
radio galaxies, iPeacock fc Nicholso find a correla- 

tion length of rp « 15.7 Mpc Cfe=0.7'). 'Magliocchctt i et alJ 
1,2004) calculate rp = 13.0±0.9 Mpc for 2dFGRS/FIRST 
radio galaxies with AGN signatures in their optical spectra 
(before correcting for redshift space distortions). These ra- 
dio galaxies therefore cluster with a strength between that 
of normal galaxies and rich clusters of galaxies. This is con- 
sistent with their preferential location in poor groups of 
galaxies. By deprojecti ng the angula r corre lation function 
of NVSS radio sources, iBlake fc Walll ^2002^ find a correla- 
tion length of rp ~ 8.6 Mpc. In this case, the clustering signal 
comes from radio galaxies with a median redshift of 2: ~1 
but the correlation function has been corrected for growth 
of clustering under linear theory JPeeblejIlQS'ol) . 

With the advent of large ar ea radio survey s down 
to fa inter flux densities iCondon e t al. 199^; .Becker et alJ 
I1995I) , we can deflne lower radio flux density limited surveys 
over large areas of the sky. We have defined a new radio 
galaxy redshift survey: the Texas-Oxford NVSS Structure 
(TONS) survey. This survey comprises the largest complete 
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spectroscopic sample of low luminosity (predominantly FRI) 
radio galaxies over large contiguous areas yet compiled. By 
going to fainter radio flux densities (si.4 ^3 mjy) and im- 
posmg optical limits {E i?<19.5), the TONS survey is 
optimised for looking at clustering of objects in a region of 
moderate (0^2^ 0.5) redshifts. The radio flux density 
limit of si.4 >3 mJy increases the space density of radio 
galaxies powered by AGN but does not go so deep that the 
sample becomes dominated by galaxies whose radio emission 
originates from star formation mechanisms. 

In this paper, we determine the clustering strength of 
the low luminosity TONS radio galaxies at moderate red- 
shifts. From the correlation length we will be able to esti- 
mate the bias of this population and compare it to previous 
estimates of both higher luminosity radio galaxies and sim- 
ilar radio galaxies in the local Universe. By cross-matching 
radio surveys with large spectr oscopic surveys such a s the 
Sloan Digital Sky Survey fS pSS:lAbazaiian et allbOO^) and 
2dF galaxy redshift survey JColless et aLlbOOlT) with NVSS 
and FIRST, even larger samples can now be compiled (see 
e.g., Jarvis, Clewley & Brand in prep.). 

The paper will be structured in the following way: In 
Sec. 2 we summarise the selection techniques used in the 
Texas-Oxford NVSS Structure (TONS) redshift survey of 
radio galaxies: this technique has al ready been discuss ed in 
detail for the TONS08 sub-sample (Brand et al.'200^. We 
present the full TONS12, TONSOSw and T0NS16w sub- 
samples and their spectra in Sec. 3. Sec. 4 summarises the 
method used to obtain a model redshift distribution and 
compares this to the observed redshift distribution of each of 
the TONS sub-samples. In Sec. 5 and Sec. 6, we describe the 
method used to calculate the two point correlation function 
and present the results for the different sub-samples. Sec. 7 
is a discussion which compares the results to previous work. 
The main conclusions are summarized in Sec. 8. 

Unless otherwise stated, we assume a spatially flat 
ACDM Universe throughout the paper with the following 
values for the cosmological parameters: Hubble constant: 
Ho = 70 km s"^Mpc"^; h = ifo/lOO = 0.7; matter density 
parameter at z=0: r2M(0) = 0.3; vacuum density parameter 
at z=0: D.AiO) = 0.7. 

2 SAMPLE SELECTION 

The Texas-Oxford NVSS Structure (TONS) survey com- 
prises three independen t regions on the sk y selected in the 
same areas as the 7CRS ("Willott et al. 200^and the TexOx- 
1000 (TOOT) survey LHiU fc Rawlines 200a) . The TONS08 
(Texas -Oxford NVS S Struc ture 08'') sub-sample is the sub- 
ject of 'Bran d etal] (|20o3). The TONS12, TONS16w and 
TONS08w surveys are presented in this paper. The sky re- 
gions surveyed by these sub-samples are shown in Table. 
Table. |5| summarises the properties of the different TONS 
sub-samples in addition to the various radio galaxy sam- 
ples that are used or referred to throughout this paper. 
TONS08 and TONS12 have the faintest radio flux densities 
limit (si.4 ^3 mJy), while TONS08w and T0NS16w were 
designed to survey larger areas but due to limited telescope 
time, were less radio deep (si.4 >30 mJy). TONS08w was 
specifically designed to determine the spatial extent of the 



TONS08 super-structur es which appeared to be unbounded 
by the TONS08 survey jBrand et alJl2003^ . 

Unlike the low-frequency selected 7CRS and TOOT, 
the TONS survey is selected at 1.4 GHz from the NVSS. 
For objects of typical spectral index, the TONS08 and 
TONS 12 sub-samples go to fainter radio flux densities than 
TOOT (which has a 151 MHz flux density sisi limit of 100 
mJy, corresponding to S1.4 «20 mJy for radio spectral in- 
dex'^Q ~0.8). In addition, TONS has an optical magnitude 
limit {E ^ R ~19.5) and colour cut {B-R >1.8) imposed 
on it. This means that we can efficiently trace large-scale 
structure at moderate {z <0.5) redshifts using AGN-fuelled 
radio galaxies. 

Full details of the sample selection, observations and 
data reduct i on of the TONS08 survey can be found in 
iBrand et all J2003l) . TONS12, TONS16w and TONS08w 
were all selected, observed and reduced using the same pro- 
cedure. The only exception to this is in the selection of the 
TONS08w sample: sample members were selected by only 
considering candidates with an angular distance of ^ 10 arc- 
sec between the optical and radio positions. This sample 
may therefore suffer from some incompletenesses. 

The selection criteria for all TONS surveys were based 
on cross-matching positions of objects in radio and opti- 
cal s urveys. An initial se lection was made of NVSS tar- 
gets iCondon et alJ [l9^8^ in the chosen area of sky with 
1.4 GHz flux densities, S1.4 >3 mJ y. The NVSS positions 
were matched with APM positions llMcMahon et aljl20o3) 
of objects with E <19.5. We selected any objects with APM 
and NVSS positions with an offset of ^ 20 arcsec from each 
othe r. We then plo tted radio contours from the FIRST sur- 
vey (iBecker et alj ll'9951 over optical POSS-II images and 
overplotted the positions of the APM and NVSS objects, 
identifying real identifications b y eye. For full details of this 
technique we refer the reader to lBrand et al.l ll2003ll . 

To correct for plate-to-plate variations in the magni- 
tudes obtained from the APM survey, we modified the O 
and E magnitudes using corrections derived from compari- 
son of the APM magnitudes to GSC-2 magnitudes for each 
POSS-II plate (R. White priv. com.). We modified the mag- 
nitude cut from our original target of 19.5 to 19.83 and 19.64 
in the TONS08 and TONS12 sub-sample respectively (the 
magnitude corresponding to the new limit of the most shal- 
low plate). All following analysis uses these corrected values. 

The final selection criteria was that O — E J5I.8. This 
cuts out all the bluer objects which tend to be stars (which 
by coincidence lie close to the line of sight to the radio 
source), quasars (which can lie at much higher redshifts) 
and star-burst galaxies (which tend to be at low redshifts 
due to their low radio luminosity). 

Table. [^summarises the total number of objects in each 
of the TONS sub-samples as well as the predicted number 
from a model redshift distribution described in Sec. 14.11 



3 OBSERVATIONS 

Optical spectra were obtained during the period October 
2000 - May 2003 on the 2.6m Nordic Optical Telescope 

^ a is the spectral index for radio sources where the radio flux 
density oc u~°', where u is the observing frequency 
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sample RA Dec 

TONS08 08'"10'"20= ^ RA 08''29'"20= 24°10'00" < DEC ^ 29°30'00" 

TONS12 12'*48'"00= ^ RA ^ 13''12'"00= 33°00'00" < DEC ^ 38°00'00" 

TONS08W 08'*00™00'' ^ RA ^ 08''41'"00= 21°50'00" < DEC ^ 31°50'00" 

TONS16W 16'"24™00'' ^ RA < 16''56'"00'' 42°00'00" ^ DEC < 51°00'00" 



Table 1. The RA and DEC limits for the different TONS sub-samples. 



Survey Type radio limits optical limits colour cut area (sr.) ref. 



TONS08 


ROz 


si.4 >3 


E <19.83 


B 


- R >1.8 


0.00688 


iBrand et al.l 1200311 


TONS12 


ROz 


si.4 >3 


E <19.64 


B 


- R >1.8 


0.00744 


this paper 


TONS16W 


ROz 


si.4 >30 


E <19.5 


B 


- R >1.8 


0.0151 


this paper 


TONS08W 


ROz 


si.4 >30 


E <19.5 


B 


- R >1.8 


0.02786 


this paper 


Lacy 


ROz 


si.4 >20 


17< E <20.2 




APM5 


0.0122 


iLacvl l200Cri 


TOOT 


Rz 


S151 >100 


NONE 




NONE 




iHill & Rawlinffd |200,* 


7CII 


Rz 


S151 >500 


NONE 




NONE 




IWillott et alj I2002D 


Sadler 


ROz 


si,4 >2.8 


14< Bj <19.4 




AGNi 


0.099 


ISadler et al.l I2002H 


2dF 


Oz 


NONE 


Bj <19.4 




NONE 


0.61 


IColless et al.l |2001|1 


NVSS2 


R 


si.4 >2.5 


NONE 




NONE 


10.36 


ICondon et alj Il996tl 


FIRST^ 


R 


si.4 >2 


NONE 




NONE 


2.61 


iBeclser et al.l ll99Sf) 


PN 


Rz 


si.4 >500 


NONE* 




NONE 


9.3 


iPeacocli & Nicholsonl 1199111 



Table 2. Table summarising the various samples used in and/or referred to in this paper. R denotes radio galaxy surveys, O denotes 
surveys with an optical magnitude cut and z denotes spectroscopic redshift surveys. The radio limits are in mjy. Notes: 1. AGN were 
identified using Principal Component Analysis (PCA) applied to the spectra iF olkes et alll99gl) . 2. NVSS is 50 pe r cent complete at 2.5 
mJy a nd 99 per cent complete at 3.5 mJy. 3. FIRST is 95 per cent complete at the given limiting flux densities. 4. iPeacock fc Nicholsonl 
^1993) impose a redshift limit of 2=0.1. 5. Galaxies were identified from APM charts. 



sample total number predicted number 



TONS08 


84 


95.1±4.5 


TONS12 


107 


97.6±4.5 


TONSOSw 


47 


74.8±3.5 


TONS16W 


40 


40.6±1.9 


Lacv ('2000') 


36 


47.5±2.7 


Sadler et al. (2002) 


433 


428.0±27.9 


TONS0812 


191 


192.7 ±6.4 


TONS0812.SS 


92 


46.9 ±3.4 


TONS0812.nonSS 


99 


145.8±5.4 



Table 3. The total number of radio galaxies in each of the 
TONS sub-samples compared to that predicted by the model 
redshift distribution. TONS0812 is the combined TONS08 and 
T0NS12 sub-samples, TONS0812_SS and TONS0812_nonSS are 
the combined TONS08 and TONS12 sub-samples within super- 
structure regions and out of super-structure region respectively 
(see Sec.lH^. 

(NOT) using the Andalucia faint object spectrograph, the 
4.2m Wilham Herschel telescope (WHT) using ISIS, the 
2.7m Smith reflect or at McDonald with the Imaging grism 
instrument (IGI) jHrneiai] |20o3) , and the Hobby- Eberly 
telescope (HE X) using the Ma rcario low resolution spectro- 
graph (LRS) jHiU et alJll99l . 

The spectra were reduced using standard IRAF proce- 
dures. Spectra and re dshifts for the TON S08 sub-sample 
have been presented in lBrand et alJ i2003l) . Spectra for the 
observed radio galaxies (with no previously known redshift) 
for the T0NS12, TONSOSw and TONS16w sub-samples are 
shown in Fig. Fig. |5| and Fig. |3 respectively. As found 



for the TONS08 sub-sample, the vast majority of objects in 
the entire TONS survey were identified as moderate redshift 
radio galaxies. All spectra are smoothed for presentation 
purposes and sorted by redshift for clarity. In most cases, 
redshifts were determined from absorption lines (very few 
objects exhibited emission lines). The estimated redshifts 
are presented in Table. |1| Table. |K|and Table. |S| These were 
checked by performing a cross-correlation between the spec- 
tra and both a rest frame composite sp ectrum obtained from 
the TONS08 data jBrand et alj|2003h and an ev olved pop- 
ulation GISSEL model (iBruzual fc Chariolll993^ . 

Of all the sub-samples, only TONSlGw overlaps with 
the sec ond data release of the Sloan Digital Sky Survey 
(SDSS: lAbazaiian et aljl2 bo4l. 11/40 of the T0NS16w sam- 
ple have spectroscopic redshifts in SDSS. In all cases, the 
redshifts agree to within 0.003. 



4 THE TONS REDSHIFT DISTRIBUTIONS 
4.1 Obtaining a universal N(z) 

To perform clustering analysis on a survey, it is crucial to 
know the underlying redshift distribution of the population 
under study. This is often obtained by simply using a red- 
shift distributio n that is either sm oothed or fitted to the 
data itself (e.g. ISteidel et aLllToQ^ . However, the obvious 
presence of large-scale structure in our sub-samples means 
that this will not be representative of the true underlying 
redshift distribution. 

We required a model redshift distribution derived from 
a data set over a sufficiently large area to reduce the effects 



a. 



O CO 

a £L 

a. ^ 2 

° 2 H 

^ CT- 

P o ^ 

B' p 

S <T> 2- 



a 



a ii^ 

m B • 

H & > 

o S s 



03 



S. 0; ^ 



a 

O 

a^ ^ 

ft) Q ft) 
<D 03 <t) 



^ > 

03 ^ 

fl^ <t) 

w 0) 



O TO 



fi, 0) 



a" 

o m 

O ^ 

m (Z) 



■0 ■ 



^ a. 
e- 2 



TO 



0) 



o 
3 



5 3 

, p g g 

03 O ^5 

M ft) ^ 

. ^ » o 

is §- ^ 

d 03 

03 2 

pj O P 

S a o 

g » P 

^ ^ 

CD C/: g 

£. 5 ^ 

S s ?= 

S " z 

o e o 

S &. 03 

? ^ ^ 

^ 

o n pi 

^ TO H 

- I O 

^ ^ I 

M 03 

C/3 3 ^ 

El ^ ^ 



s ^ 

p ft) 

0) £ 

CD P". 



0) 
O. ^1. 



CD 



:;■ O 
/) O 



i-i a 
P TO 



O. (D S- 

03 h-. 

□ ^ B 

3"' «2 B" 



Ki-imo 


T 


)|0SfO,,< 


Diilo 




T!A(J 


2000) 


r)E("'(.72000) 






F 









/■) 




Kolop 


TONS12_03Q 


2 


7 


2QQ3 


May 


12 


48 


32.94 


34 


19 


01.85 


10.1 


7 


76 


17.96 


21.05 


0.234 





2317 


0.2313 




TONS12_045 










12 


48 


41.37 


33 


00 


55.50 


40.2 


38 


23 


19.46 


22.05 


1.23 








TOOT12_100:105 


TONS12_048 


2 


7 


2QQ3 


May 


12 


48 


45.58 


34 


46 


38.50 


3.5 






18.49 


21.56 


0.375 





3715 


0.3747 




TONS12.063 










12 


49 


01.49 


36 


15 


35.39 


10.2 


3 


79 


19.25 


22.05 


0.560 








TOOT12.100:112 


TONS12_112 


2 


7 


2QQ3 


May 


12 


49 


58.72 


36 


21 


25.16 


16.1 


5 


38 


18.74 


21.76 


0.397 





3954 


0.3968 




TONS12.116 


2 


7 


2QQ3 


May 


12 


50 


04.67 


36 


39 


33.22 


5.4 


3 


99 


16.04 


19.13 


0.161 





1620 


0.1613 




TONS12_122 










12 


50 


08.65 


33 


09 


37.83 


7.4 


3 


01 


6.65 


10.71 


0.024 










TONS12.130 


2 


7 


2003 


May 


12 


50 


16.76 


34 


45 


16.30 


5.8 


4 


41 


17.13 


20.11 


0.338 





3378 


0.3384 




TONS12_134 


2 


7 


2003 


May 


12 


50 


24.98 


35 


55 


08.97 


6.3 


4 


10 


18.30 


21.70 


0.323 





3244 


0.3258 




TONS12.145 


2 


7 


2003 


May 


12 


50 


37.13 


35 


00 


08.69 


30.6 


31 


71 


19.00 


21.72 


0.450 





4492 


0.4522 




TONS12_150 


2 


7 


2003 


May 


12 


50 


44.14 


33 


43 


22.19 


3.7 


4 


13 


18.32 


21.49 


0.220 





2203 


0.2170 




TONS12_151 










12 


50 


44.94 


35 


24 


00.69 


3.6 


6 


87 


12.56 


15.46 


0.032 








2MAS 


TONS12_165 


2 


7 


2003 


May 


12 


50 


57.68 


35 


46 


52.29 


14.4 


8 


56 


17.99 


20.53 


0.321 





3218 


0.3222 




TONS12_189 


2 


7 


2003 


May 


12 


51 


16.87 


36 


11 


37.23 


5.7 






18.57 


20.85 


0.318 





3183 


0.3196 




TONS12.211 


2 


7 


2003 


May 


12 


51 


42.33 


37 


29 


53.59 


9.4 


6 


21 


17.99 


21.48 


0.307 





3074 


0.3078 




TONS12_218 


2 


7 


2003 


May 


12 


51 


47.08 


37 


12 


04.56 


4.3 


2 


63 


19.31 


21.48 


0.420 





4199 


0.4185 




TONS12.226 


2 


7 


2003 


May 


12 


51 


53.46 


35 


19 


24.58 


15.2 


3 


04 


17.17 


20.58 


0.407 





4071 


0.4070 




TONS12_228 


2 


7 


2003 


May 


12 


51 


53.90 


37 


27 


43.04 


30.6 


17 


62 


17.91 


20.92 


0.290 





2896 


0.2917 




TONS12.252 


2 


7 


2003 


May 


12 


52 


29.72 


35 


13 


44.49 


3.8 


2 


42 


18.33 


21.62 


0.407 





4074 


0.4068 




TONS12_26Q 










12 


52 


40.25 


33 


10 


58.43 


263.3 


222 


66 


19.57 


22.05 


0.509 








TOOT12_300:178 


TONS12.27Q 


2 


7 


2003 


May 


12 


52 


48.45 


37 


42 


06.78 


11.2 


7 


97 


18.87 


21.25 


0.353 





3530 


0.3531 




TONS12_271 


2 


7 


2003 


May 


12 


52 


48.51 


33 


36 


21.24 


19.0 


17 


33 


18.93 


22.05 


0.460 





4587 


0.4612 




TONS12_272 


2 


7 


2003 


May 


12 


52 


54.18 


35 


37 


20.53 


5.8 


4 


36 


17.90 


21.65 


0.322 





3211 


0.3206 




TONS12.277 


2 


7 


2003 


May 


12 


52 


58.35 


35 


12 


45.87 


12.5 


13 


52 


18.59 


22.05 


0.247 





2470 


0.2479 




TONS12_286 










12 


53 


07.69 


33 


29 


20.31 


3.8 


2 


27 


12.04 


17.01 


0.129 








2MASS 


TONS12.287 


2 


7 


2003 


May 


12 


53 


07.95 


33 


46 


05.39 


21.3 


20 


76 


18.96 


22.05 


0.422 





4225 


0.4226 




TONS12_301 


2 


7 


2003 


May 


12 


53 


24.30 


34 


35 


12.32 


362.0 


352 


94 


18.67 


22.05 


0.034 


- 




- 


TOOT12_500:192 


TONS12.306 


2 


7 


2003 


May 


12 


53 


33.03 


36 


43 


39.71 


3.1 


2 


62 


17.06 


21.01 


0.306 





3042 


0.3040 




TONS12_312 


2 


7 


2003 


May 


12 


53 


39.40 


35 


11 


06.26 


10.9 


9 


19 


19.26 


22.05 


0.458 





4565 


0.4587 




TONS12_326 


2 


7 


2003 


May 


12 


53 


59.25 


37 


37 


55.18 


3.6 


4 


11 


18.88 


21.61 


0.308 





3075 


0.3068 




TONS12_331 










12 


54 


11.91 


36 


30 


59.54 


72.1 


62 


95 


17.42 


20.86 


0.217 








Benn 


TONS12_338 


2 


7 


2001 


Apr 


12 


54 


20.30 


36 


05 


00.04 


7.3 


4 


64 


16.49 


20.03 


0.223 





2233 


0.2232 




TONS12_373 










12 


54 


59.11 


34 


00 


45.74 


4.6 


5 


78 


18.21 


21.42 


0.150 








TOOT12_300:223 


TONS12_400 


2 


7 


2002 


Feb 


12 


55 


17.25 


35 


53 


58.30 


60.5 


46 


76 


18.11 


22.05 


0.364 





3631 


0.3634 




TONS12.413 


2 


7 


2001 


Apr 


12 


55 


30.17 


35 


46 


03.65 


3.4 


3 


52 


16.49 


20.78 


0.327 





3273 


0.3278 




TONS12_419 










12 


55 


33.43 


36 


57 


52.22 


139.7 


104 


05 


19.00 


22.05 


0.206 








Benn 


TONS12.421 










12 


55 


36.87 


34 


56 


46.96 


59.0 


48 


79 


13.38 


18.19 


0.186 








Benn 


TONS12_425 


2 


7 


2001 


Apr 


12 


55 


39.42 


36 


32 


31.21 


6.3 





94 


16.66 


20.32 


0.221 





2202 


0.2200 




TONS12.427 


2 


7 


2002 


Feb 


12 


55 


45.82 


34 


31 


59.02 


4.9 






19.62 


22.05 


0.426 





4261 


0.4263 




TONS12_443 


2 


7 


2002 


May 


12 


56 


00.71 


36 


23 


47.90 


3.1 


2 


12 


15.06 


18.20 


0.102 


- 




- 




TONS12_446 


2 


7 


2001 


May 


12 


56 


09.16 


35 


38 


11.31 


16.5 


16 


73 


11.31 


14.81 


0.076 





0768 


0.0780 




TONS12_453 










12 


56 


13.09 


36 


38 


33.87 


37.3 


21 


07 


19.26 


22.05 


0.467 








TOOT12_500:243 


TONS12_469 


2 


7 


2002 


Feb 


12 


56 


26.58 


36 


00 


07.83 


5.5 


2 


72 


18.20 


21.01 


0.347 





3471 


0.3483 




TOKS12_511 


2 


7 


2001 


May 


12 


57 


07.37 


37 


25 


31.02 


4.4 






17.85 


22.05 


0.325 





3255 


0.3274 




TONS12_351 










12 


57 


57.42 


35 


40 


22.79 


61.2 


56 


64 


15.07 


19.82 


0.219 








Benn 


TONS12.571 


2 


7 


2001 


Apr 


12 


58 


12.76 


33 


35 


28.74 


3.7 


3 


29 


9.08 


12.84 


0.077 





0761 


0.0760 




TONS12_578 










12 


58 


27.42 


36 


58 


59.98 


104.5 


66 


28 


19.40 


22.05 


0.436 








TOOT12_500:289 


TONSI2.088 










12 


58 


34.68 


36 


51 


14.65 


26.3 


17 


69 


18.82 


22.05 


0.464 








Benn 


TONS12_o93 


2 


7 


2002 


Mai- 


12 


58 


42.77 


34 


03 


57.62 


25.0 


21 


41 


16.23 


20.14 


0.248 





2471 


0.2469 




TONS12.597 










12 


58 


51.64 


35 


02 


44.31 


26.6 


34 


05 


17.23 


20.99 


0.248 








Benn 


TONS12_601 


2 


7 


2001 


May 


12 


58 


56.26 


37 


23 


56.57 


6.5 


8 


09 


8.43 


11.60 


0.037 





0371 


0.0365 




TONS12.6Q7 










12 


59 


00.26 


34 


50 


48.14 


12.3 


6 


01 


8.36 


22.05 


0.003 








Benn 


TONS12_618 










12 


59 


08.92 


37 


18 


36.64 


26.5 


9 


23 


0.45 


7.77 


0.016 








Benn 


TONS12.619 


2 


7 


2001 


May 


12 


59 


09.36 


36 


33 


43.15 


6.6 


3 


30 


15.01 


19.81 


0.218 





2173 


0.2174 





Co 

a 

s 



Name 



Tolof.™|)0 Dale 



T{A(.72000) r)E("'(.72000) ,s'n\-ss 



TONS12_638 


2.7 


2002 


Feb 


12 


59 


35, 


.48 


37 


13 


42, 


.39 


5, 


.3 


TONS12.692 


2.7 


2002 


Feb 


13 


00 


25, 


.01 


33 


34 


30, 


.7 


7, 


.2 


TONS12_726 










13 


01 


06, 


.41 


35 


40 


26, 


.68 


54, 


.7 


TONS12.731 


2.7 


2002 


Feb/Mar 


13 


01 


10, 


.97 


35 


54 


32, 


.63 


8, 


.7 


TONS12.745 










13 


01 


25, 


.24 


36 


58 


08, 


.55 


28, 


.7 


TONS12.761 


2.7 


2002 


Apr 


13 


01 


45, 


.53 


35 


30 


08, 


.35 


5, 


.1 


TONS12_775 


HET 


2002 


Feb 


13 


02 


00, 


.81 


37 


40 


19, 


.31 


6, 


.2 


TONS12.818 


2. 


.7 


2002 


Feb 


13 


02 


57, 


.79 


33 


29 


54, 


.77 


5, 


.7 


TONS12_827 










13 


03 


10 


.59 


34 


57 


46 


.89 


9, 


.3 


TONS12.829 


2. 


.7 


2002 


Feb 


13 


03 


12 


.22 


36 


51 


58 


.05 


5, 


.5 


TONS12_855 










13 


03 


43 


.94 


37 


56 


09 


.90 


575 


.9 


TONS12_866 


2. 


.7 


2001 


May 


13 


03 


55 


.53 


34 


30 


46 


.00 


7, 


.4 


TONS12_876 










13 


04 


07 


.29 


37 


09 


07 


.63 


319 


.6 


TONS12_895 


2, 


,7 


2001 


May 


13 


04 


26 


.77 


33 


26 


50 


.14 


4, 


.8 


TONS12.925 


2. 


.7 


2002 


Feb 


13 


05 


01 


.40 


36 


12 


51 


.62 


4, 


.2 


TONS12_958 


2. 


.7 


2002 


Feb 


13 


05 


34 


.23 


33 


39 


01 


.34 


7, 


.4 


TONS12.962 


2, 


.7 


2001 


May 


13 


05 


37, 


.83 


36 


51 


38 


.28 


8, 


.2 


TONS12_971 










13 


05 


51, 


.68 


36 


39 


28, 


.51 


133, 


.8 


TONS12.972 










13 


05 


51, 


.79 


36 


22 


53, 


.41 


208, 


.9 


TONS12.974 


2. 


.7 


2002 


Feb 


13 


05 


55, 


.62 


33 


33 


08, 


.12 


4, 


.7 


TONS12.992 


2. 


.7 


2001 


May 


13 


06 


06, 


.65 


37 


49 


15, 


.56 


4, 


.1 


TONS12_993 










13 


06 


07, 


.19 


33 


38 


13, 


.22 


20, 


.7 


TONS12_994 


2. 


.7 


2002 


Feb 


13 


06 


08, 


.26 


34 


10 


39, 


.54 


19, 


.5 


TONS12.999 










13 


06 


14, 


.13 


35 


30 


06, 


.09 


66, 


.9 


TONS12_1004 










13 


06 


18, 


.24 


34 


22 


22, 


.78 


7, 


.1 


TONS12.1018 










13 


06 


24 


.24 


36 


42 


19 


.46 


17 


.5 


TOKS12_1035 










13 


06 


40 


.04 


34 


22 


19 


.94 


16 


.2 


TONS12.1041 


2, 


,7 


2002 


Feb 


13 


06 


45, 


.74 


33 


07 


47, 


.28 


7 


.7 


TONS12_1057 


2, 


,7 


2002 


Feb 


13 


07 


10 


.61 


36 


56 


26 


.98 


14, 


.3 


TONS12.1061 


2. 


.7 


2001 


May 


13 


07 


11 


.58 


34 


54 


08 


.23 


3 


.5 


TOKS12_1062 










13 


07 


12 


.50 


34 


39 


38 


.15 


10 


.3 


TONS12_1068 


2, 


,7 


2002 


Feb 


13 


07 


23 


.22 


36 


14 


30 


.93 


6 


.0 


TONS12_1080 


2. 


.7 


2002 


Feb 


13 


07 


42 


.35 


37 


44 


19 


.32 


3 


.8 


TONS12_1099 


2. 


.7 


2002 


Mar/Apr 


13 


08 


01 


.00 


36 


22 


48 


.53 


5 


.0 


TONS12.1112 










13 


08 


17, 


.69 


33 


23 


55, 


.94 


50, 


.5 


TONS12_1139 


2. 


.7 


2001 


May 


13 


08 


54, 


.89 


34 


15 


49, 


.74 


3, 


.6 


TONS12.1155 


2. 


.7 


2001 


May 


13 


09 


15, 


.91 


36 


48 


56, 


.80 


12, 


.3 


TONS12_118Q 


2. 


.7 


2001 


May 


13 


09 


40, 


.13 


34 


11 


33 


.61 


3, 


.6 


TONS12.1196 










13 


09 


53, 


.89 


33 


59 


28, 


.05 


68, 


.3 


T0NS12_12Q1 










13 


10 


00, 


.03 


33 


56 


59, 


.10 


29, 


.9 


TONS12.1207 


2. 


.7 


2003 


May 


13 


10 


05, 


.10 


37 


28 


57, 


.19 


15, 


.2 


TONS12_1211 


2. 


.7 


2003 


May 


13 


10 


08, 


.06 


35 


10 


18, 


.92 


13, 


.2 


TONS12_1231 


2. 


.7 


2003 


May 


13 


10 


30 


.55 


33 


20 


16 


.69 


26 


.0 


TONS12.1242 


2. 


.7 


2003 


May 


13 


10 


39 


.50 


34 


55 


44, 


.94 


4, 


.8 


TONS12_1245 


2. 


.7 


2003 


May 


13 


10 


44, 


.75 


34 


55 


25 


.18 


10 


.2 


TONS12.1255 


2, 


,7 


2003 


May 


13 


10 


51 


.49 


35 


10 


31 


.02 


10 


.5 


TONS12_1267 


2. 


.7 


2003 


May 


13 


11 


00 


.15 


36 


38 


19 


.21 


3 


.9 


TONS12.1283 










13 


11 


22 


.41 


34 


25 


02 


.13 


22 


.8 


TONS12_1292 










13 


11 


31 


.16 


36 


16 


54, 


.94 


11 


.4 


TONS12.1298 


2. 


.7 


2003 


May 


13 


11 


38 


.45 


35 


51 


55 


.91 


13 


.6 


TONS12_1306 










13 


11 


46 


.01 


34 


33 


18 


.40 


11 


.9 


TONS12.1311 


2. 


.7 


2003 


May 


13 


11 


50, 


.93 


34 


56 


32, 


.79 


5, 


.8 


TONS12_1313 


2. 


.7 


2003 


May 


13 


11 


52, 


.87 


34 


59 


19, 


.83 


21, 


.0 





F 


(") 




( OPi ) 






Kol op 




5.71 


18.35 


22.05 





.369 


0.3693 


0.3694 






4.83 


15.73 


18.17 





.094 










53.04 


19.25 


22.19 





.622 






TOOT12_100:353 




8.05 


18.51 


22.19 





.456 


0.4560 


0.4568 






16.07 


18.36 


22.19 





.424 






Benn 




2.92 


18.16 


20.61 





.437 


0.4368 


0.4353 






5.89 18.79 


22.19 







.328 


0.3266 


0.3266 






7.26 


18.79 


21.87 





.390 


0.3893 


0.3911 






9.54 


17.40 


20.23 





.565 






Benn 




5.41 


18.26 


22.19 





.423 


0.4223 


0.4229 






569.48 


18.31 


20.30 





.470 






Benn 




3.02 


17.95 


21.30 





.283 


0.2840 


0.2839 






304.89 


19.49 


21.09 





.425 






Benn 




3.98 


16.47 


19.71 





.200 


0.1993 


0.1998 






4.38 


18.51 


21.80 





.373 


0.3716 


0.3727 






2.29 


17.53 


20.53 





.242 


0.2422 


0.2414 






5.67 


17.80 


21.28 





.275 


0.2750 


0.2760 




Th 


123.93 


17.31 


21.61 





.312 






Benn 


198.19 


17.61 


20.92 





.315 






Benn 


as 


4.00 18.98 


22.19 







.401 


0.4006 


0.4026 






3.27 


14.99 


16.82 





.109 


0.1082 


0.1104 






18.45 


18.15 


22.01 





.302 






Benn 


12.63 


18.42 


21.65 





.305 


0.3033 


0.3033 




dus 


64.38 


17.76 


21.93 





.373 






Benn 


5.42 


17.26 


20.91 





.301 






Benn 


tterin 


9.85 


18.02 


22.19 





.312 






Benn 


8.05 


15.67 


18.88 





.239 






Benn 


7.65 


18.71 


22.19 





.302 


0.3023 


0.3052 






6.46 


18.16 


22.19 





.366 


0.3660 


0.3665 






3.99 


16.85 


20.04 





.242 


0.2414 


0.2424 




f ra 


6.94 


16.67 


21.58 





.317 


0.3191 


0.3200 




2.95 


18.63 


22.19 





.320 






Benn 


dio 


2.42 


18.42 


21.75 





.394 


0.3936 


0.3944 




4.19 


19.34 


22.19 





.335 


0.3346 


0.3355 






43.61 


16.14 


19.61 





.254 






Benn 




3.90 


17.49 


20.92 





.267 


0.2666 


0.2685 






7.48 


17.27 


20.89 





.239 


0.2392 


0.2386 




3.33 


17.33 


20.73 





.246 


0.2457 


0.2447 




a 


67.67 


18.58 


22.19 





.464 






Benn 




24.85 


16.96 


20.43 





.44 






TOOT12_100:552 




13.39 


14.81 


17.66 





.147 


0.1471 


0.1467 




as 


10.03 


15.22 


19.08 





.200 


0.1993 


0.1975 




26.80 


18.76 


21.70 





.434 


0.4332 


0.4306 




TONS , 


1.03 


18.38 


22.19 





.394 


0.3941 


0.3941 




8.80 


17.01 


20.73 





.287 


0.2871 


0.2873 




8.19 


18.49 


21.25 





.270 


0.2674 


0.2621 




2.04 


19.57 


22.19 





.318 


0.3190 


0.3224 




11.94 


9.81 


12.37 





.037 










6.91 


7.13 


9.09 





.004 










3.63 


15.73 


19.86 





.184 


0.1843 


0.1843 






6.76 


18.63 


21.57 





.24 






TOOT12_100:599 




5.36 


18.17 


21.47 





.303 


0.3030 


0.3041 






15.90 


17.16 


20.14 





.302 


0.3013 


0.3015 
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Name 


Tolosro]>o 


Dill 




T{.A(J2000) 










K 


( ) 




(OPi ) 








TONSQ8w_0232 


HET 


2000 


Nov 


08 


01 


44.88 


28 


35 


02 


.33 


172, 


.5 


175.15 


16.21 


20.13 





.267 


0.2678 


0.2661 




TONSQ8W.0401 


2.7 


2001 


Feb 


08 


02 


52.96 


25 


52 


55 


.25 


30, 


.3 


30.61 


12.10 


14.87 





.082 








TONSQ8W.0513 


HET 


2000 


Nov 


08 


03 


38.71 


29 


06 


57, 


.32 


165, 


.6 


152.39 


17.92 


21.43 





.302 


0.3026 


0.3017 




TONSQ8W.06Q9 


HET 


2001 


Feb 


08 


04 


22.49 


22 


16 


42, 


.36 


33, 


.8 


33.73 


18.85 


22.03 





.208 


0.2096 


0.3151 




TONSQ8W.0785 


HET 


2001 


Jan 


08 


05 


33.61 


30 


22 


45, 


.52 


59 


.8 


60.81 


18.34 


20.95 





.315 


0.3161 


0.3151 




TONSQ8W.0973 


HET 


2001 


Feb 


08 


06 


48.82 


24 


40 


21 


.16 


70, 


.8 


60.65 


18.33 


22.03 





.356 


0.3562 


0.3570 




TONSQ8w_0987 


HET 


2001 


Jan 


08 


06 


54.26 


26 


31 


17 


.22 


48, 


.4 


44.53 


19.33 


22.03 





.410 


0.4110 


0.4108 




TONS08W.1039 


2.7 


2001 


Feb 


08 


07 


17.64 


27 


47 


07 


.44 


86, 


.0 


84.72 


18.68 


21.27 





.357 








TONS08W.1062 


HET 


2001 


Jan 


08 


07 


29.66 


29 


05 


57 


.00 


99 


.7 


92.97 


15.12 


18.55 





.191 


0.1914 


0.1881 




TONS08W.1122 


WHT 


2001 


May 


08 


07 


52.36 


31 


48 


40 


,84 


110 


,6 


114.82 


16.63 


22.03 





.110 


0.1100 


0.1097 




TONS08w_1174 


2.7 


2002 


Feb 


08 


08 


17.91 


30 


46 


43 


.01 


330 


.3 


297.07 


18.37 


21.40 





.317 


0.3173 


0.3213 




TONS08W.1211 


HET 


2001 


Mar 


08 


08 


36.78 


26 


46 


36 


.55 


452 


.3 


457.07 


19.17 


22.03 





.544 


0.5449 


0.5466 




TONS08W.1242 


2.7 


2002 


Feb 


08 


08 


50.78 


27 


57 


58 


.19 


50 


.9 


47.96 


19.36 


22.03 





.424 


0.4239 


0.4238 




TONS08w_1561 








08 


10 


58.11 


28 


46 


18 


.94 


40 


.8 


27.28 


18.10 


22.03 





.267 






TOOT08_100:004 


TONS08w_1803 








08 


12 


43.90 


25 


02 


43 


.40 


55 


.7 


42.43 


17.66 


22.03 





.242 






TOOT08.200:015 


TONS08W.1951 








08 


13 


43.51 


26 


55 


09 


.41 


74, 


.2 


63.54 


17.76 


22.03 





.325 






TOOT08.200:026 


TONS08W.2137 








08 


14 


53.40 


28 


29 


40, 


.51 


149, 


.1 


78.57 


18.63 


21.46 





.256 






7C0811+2838 


TONSQ8w_2235 


WHT 


2001 


May 


08 


15 


36.51 


29 


50 


23 


.36 


62, 


.3 


51.08 


17.95 


20.29 





.336 


0.3358 


0.3339 




TONS08W.2415 








08 


16 


51.65 


28 


22 


24 


.30 


41, 


.7 


38.32 


18.99 


22.03 





.351 






TOOT08.200:049 


TONSQ8W.2445 


2.7 


2002 


Feb 


08 


17 


03.63 


29 


51 


47, 


.01 


91, 


.5 


89.58 


19.36 


22.03 





.370 


0.3718 


0.3703 




TONSQ8W.2467 


HET 


2001 


Feb 


08 


17 


11.57 
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Figure 1. Spectra of all the radio galaxies in TONS12 that are 
not in TOOT, reported by Benn (priv. comm.) and/or do not 
have redshifts from the literature. See Table. |^ for details. The 
Flux density is in units of 10"^" Wm~^A . Shown here is an 
example spectrum. Plots of all the spectra can be found in the 
version to be published in MNRAS in early 2005. 





Figure 2. Spectra of all the radio galaxies in TONSOSw that 
are not in TONS08, TOOT08 and/or not not have redshifts from 
the literature. See Table. 151 for details. The flux density is in units 
of 10"^" Wm~^A . Shown here is an example spectrum. Plots 
of all the spectra can be found in the version to be published in 
MNRAS in early 2005. 



of large-scale structure and which we could modify for differ- 
ent optical magnitude and radio flux density limits of each 
of our sub-samples. There are currently no model redshift 
distributions derived from surveys goin g down to NVSS / 
FIRST flu x densities in the li terature. iDunlop fc Peacodd 
and lWillott et all teOOlT) model the radio luminosity 



J199' 

function of relatively bright (si.4 >100 mjy) radio galaxies 
from which a redshift distribution can be extrapolated to 
lower flux densities and the correct radio frequency. However 




Figure 3. Spectra of all the radio galaxies in TONS16 that 
are not in TOOT16 and/or do not have redshifts from the lit- 
erature. See Table. |H] for details. The flux density is in units of 
10 Wm~^A . Shown here is an example spectrum. Plots 
of all the spectra can be found in the version to be published in 
MNRAS in early 2005. 



these models can display unphysical spikes at low redshifts 
and the extrapolation results in large uncertainties. These 
models were not accurate enough for our purposes. 

Instead, we use d a maximum likelihood technique 
jMarshall et al.lll983h to fit our own bivariate (radio and 
optical) luminosity function (BLF) to radio ga laxies in the 
2dF galaxy redshift survey iSadler et al.ll2002) . These data 
are over a sufficiently large area to smooth over any large- 
scale structure. The BLF can be integrated over any optical 
apparent magnitude and radio flux density limits to obtain 
a model redshift distribution for each of the TONS sub- 
samples. Brand et al. (2003) described the form of the bi- 
variate luminosity function that was modelled and the best 
fit parameters of this model. Note that because to TONS 
survey extends out only to moder ate redshift, the redshift 
evolution is not well constrained. IClewlev fc Jarvi^ (|2004) 
found that the lower luminosity (FR I dominated) AGN 
population is consistent with a constant co-moving space 
density with redshift. Fixing the BLF to a model with no 
evolution with redshift, resulted in no significant difference 
in the model redshift distribution. 

Because the TONS survey optical apparent magnitude lim- 
its were defined in the E{~ R) band but the BLF is defined 
for the bj band, we converted the apparent J?-band magni- 
tude limits to absolute R magnitude for each redshift and 
applied a colour correction of B — i?=1.52 (calculated from 
the difference in flux in the rest- frame of a template SEP 
of an evolved stellar population jBruzual fc CharlotjlTooS) 
measured in the 2 different b ands) . We assumed a colour 



correction of Bi — B 



-0.11 



The redshift distribution of the lSadler et al.l (120021) ra- 
dio galaxies is shown is Fig. |1| Overplotted is the fltted 
model redshift distribution. We note that although this sur- 
vey is is over 13 times the area of TONS08, it still exhibits 
features due to large-scale structure and this may bias the 
resulting fit of the BLF. We minimised this effect by effec- 
tively smoothing the redshift distribution (by modifying the 
redshifts of each radio galaxy by an amount drawn from a 
Gaussian distribution of FWHM 0.05). 

4.2 Errors 

We found the errors on the redshift distribution using the 
following method to allow for the deviations in the derived 
model parameters. We split the parameter space into dis- 
crete grid spaces. For each grid space, we calculated the 
maximum likelihood, S of this parameter combination and 
converted it to a probability (P cx exp( — S/2)). We then nor- 
malised the probabilities so that the total probability in all 
grid spaces added up to one and converted to a cumulative 
probability. We calculated the number of radio galaxies ex- 
pected for each redshift interval given the survey limits and 
the parameter values in that box. We then weighted this by 
multiplying by the probability of the parameter combina- 
tion. By adding up the weighted number of galaxies for each 
redshift interval, we should have obtained a very similar red- 
shift distribution to that obtained by the best fit parameter 
values. The advantage of this procedure is that if there is 
more than one parameter combination which fits the BLF 
well, these will also have an influence on the final redshift 
distribution. 

To calculate the errors, we performed a Monte-Carlo 
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Figure 4. The redshift distribution of the ISadler et all i2002 
sample with the model redshift distribution overplotted (solid 
line). The ±la errors on the model are overplotted (dashed lines). 
We plot the distribution using 16 and 55 redshift bins between 
2=0 and 2=0.5 for the top and bottom figures respectively. Even 
with 433 radio galaxies over a larger area, the redshift distribution 
is not smooth. 



simulation in which we chose 5000 random numbers between 
and 1. Because the cumulative probabilities for all param- 
eter space grid spaces were normalised to run from to 1, 
each random number could be associated with a unique pa- 
rameter combination. Those parameter combinations which 
have higher relative probabilities occupy a larger proportion 
of this probability space and will therefore have had a greater 
chance of being selected. The redshift distribution was then 
calculated for all parameter combinations that were selected. 
The minimum and maximum ±lcr errors at each redshift in- 
terval were calculated by determining where 67 per cent of 
the measurements fell. 

It is likely that other quantities will also have con- 
tributed to the error in the measurement of the redshift 
distribution. Errors may have been introduced by incom- 
pletenesses in both the radio and optical surveys at faint 
flux densities: The NVSS catalogue is ~ 90 per cent com- 
plete at si.4=3 mjy iCondon et al.lll998ll . Before correcting 



for plate-to-plate variations, the A PM magnitudes have a 
global rms uncertainty of 0.5 mag jMcMahon et al]l2002l) . 
We assumed an error of 0.3 mag, which is the difference be- 
tween the magnitude corrections of the two POSS-II plates 
in the TONS08 region. We incorporated this into the to- 
tal error on each redshift bin by determining the R mag- 
nitude limit in the Monte Carlo simulations as a Gaussian 
distributed value with a standard deviation of 0.3 about the 
magnitude limit. 



4.3 A comparison of the TONS sub-samples with 
their model redshift distributions 

Table. |H] shows the number of radio galaxies observed and 
predicted by our model redshift distribution in all the TONS 
sub-samples. These numbers are significantly different in 
most cases. However, it is not surprising that there are large 
variations in these numbers due to the presence of super- 
structures and voids which occur on similar scales to that of 
the survey size (~100 Mpc diameter at z ~0.3). 

Fig.|5]shows the redshift distributions overplotted with 
the model redshift distributions for the TONS08, TONS12, 
TONSOSw and TONS16w sub-samples. In aU cases, large- 
scale structure has changed the redshift distribution of our 
sub-sample significantly from that expected for a typical 
area of sky. 

Within the TONS08 sub-sample, there are two promi- 
nent spikes at z f»0.27 and z ~0.35 with, perhaps, void- 
like regions at lower redshifts, accounting for the low overall 
number of radio galaxies. Analysis of these redshift spikes; 
their significance and their implications for structure for- 
mation theories and/or radio galaxy bias is performed in 
[Brand ct ah (_2003). The TONS 12 survey covers a similar 
area to that of TONS08 but in a different region of sky. 
It too shows two redshift peaks above that expected from 
the model redshift distribution (but at redshifts z «0.23 and 
z ~0.31). That redshift spikes are found in both sub-samples 
is important as it demonstrates that redshift spikes are prob- 
ably a universal feature of radio galaxy distributions. 

The TONSOSw and 16w surveys are wider area sam- 
ples with higher flux density limits. Within the TONSOSw 
sub-sample, there are fewer radio galaxies than we would 
expect. However, this may simply be due to incomplete- 
nesses in this survey. We may have missed sources with a 
large offset between the radio and optical positions - some- 
thing more likely for higher power radio galaxies. Within 
the T0NS16w sub-sample, we observed a similar number of 
radio galaxies to that predicted from the model redshift dis- 
tribution. Although there is an excess of radio galaxies at a 
redshift 2=0.2-0.3, the number density of radio galaxies is 
insufficient to determine if there is a significant overdensity. 
This is perhaps not surprising given the lack of evidence for 
significant redshift spi kes in previous sur veys with a similar 
fiux density limit (see lBrand et alJl2003l sec. 5.3). 



5 RADIO GALAXY CLUSTERING 

5.1 Calculating the two point correlation function 

The spatial two-point correlation function is the traditional 
tool for measuring the strength of clustering and hence pro- 
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Figure 5. The redshift distribution of the TONS08, T0NS12, TONSOSw and T0NS16w sub-samples (top left to bottom right) with the 
model redshift distribution ovcrplotted (solid line). The ztlcr errors on the model are also shown (dashed lines). We plot the distribution 
using 28 redshift bins between 2=0 and 2=0.5. 



vides an easy comparison with other surveys. It is also per- 
haps the most intuitive and simple measure of clustering. It 
provides fundamental information about the galaxy distri- 
bution in the sense that it is the Fourier transform of the 
power spectrum of the density fluctuations and (when the 
clustering is still in the linear regime and the bias is known) 
can be directly related back to fluctuations in the primordial 
density fleld. The two-point correlation function measures 
the excess probability over random AP of finding a pair of 
obje cts in two volumes dVi an d AV2 with a separation, r 
(e.g. iPeacock fc NicholsonllTool : 



dP = pg[l-he(r)]dVidV2, 



(1) 



where po is the average number density. The form of the cor- 
relation function is generally well approximated by a power- 
law: 



(2) 



where ro is the correlation scale length (the separation at 
which ^=1) and 7 is the slope of the correl ation function and 
is usu ally assumed to be equal to 1.8 (e.g. lGroth fc Peebles! 
Il977tl . Thus, the strength of any clustering can be estimated 



by obtaining a value for ro. The larger the value of ro, the 
greater the clustering signal of the sample. 

In practise, the correlation function is calculated using 
an estimator. There are several different forms of this esti- 
mator, but they all involve counting pairs of objects within 
the real catalogue, a random catalogue and between the two 
at different separations, r. Korschcr ct al. (2000) find that on 
small (r~4.4/i~^Mpc) scales, the estimators are comparable. 
However, for larger sca les (r«115/i~^Mpc), the estimator of 
iLandv fc Szalav^^l993^ significantly outperforms the rest: 



DDjr) - 2DR{r) + RR{r) 
RR{r) 



(3) 



where DD(r) is the number of data-data pairs, RR is the 
number of random-random pairs and DR is the number of 
data-random pairs all calculated for separation r. This is the 
estimator that we have used for this analysis. 

In order to calculate the two-point correlation function, 
a set of random catalogues was constructed. These cata- 
logues allow a comparison of the data to randomly dis- 
tributed objects on the sky. The random catalogues were 
all constructed with the same number of objects as the real 
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catalogue and the same selection function. They were all as- 
signed a random RA and DEC within the sample boundaries 
and a random redshift selected from the redshift distribution 
expected for objects in the survey. We used the redshift dis- 
tribution calculated in Sec. 14.11 The RR and DR pairs used 
in the calculation of the correlation function were obtained 
by averaging 1000 simulations of the random catalogue to 
reduce to negligible any errors due to statistical variations 
within the random samples. 



5.2 Corrections to the two-point correlation 
function 

Small survey volumes that are unrepresentative of the over- 
all field must be corrected for the "integral constraint" . Be- 
cause the TONS survey sizes are comparable to scales at 
which the number density of galaxies in the volume can sig- 
nificantly differ from field to field, we need to account for 
the fact that they are not representative of the general field. 
Because the correlation function is only calculated for a fi- 
nite survey volume, in practise the observed number density 
n is used instead of the global number density < n >. 

Unlike most studies, we have corrected for the sample 
being under- or over-dense due to large-scale structure. This 
is possible because, in calculating the model redshift distri- 
butions for each subsample, we have already estimated the 
underlying number density of radio galaxies in the TONS 
sub-samples. The observed correlation function, ^(r)obs was 
modified to the true correlation function, ^(r)truc by a nor- 
malisation term, n/ < n >: 



(l + CMt.-uo) = (l + e(r)obs) 



(4) 



JPeacocklll999bl) where n is the total number of objects in 
the sample and < n > is the total number of objects ex- 
pected in the sample as calculated from our model redshift 
distribution (see Table. I^J. In an under-dense region, the 
number of random-random pairs calculated in the correla- 
tion function will be under-estimated and hence the correla- 
tion function will be over-estimated and must be corrected 
downwards. Applying this formula is essentially the same 
as changing the number of random objects in equation. |3 
from the number of objects in the sample to the number of 
objects expected. 

The observed correlation function is also modified be- 
cause it is measured not in real space, but in redshift space. 
In practise, the redshift measurements will be affected by 
large-scale coherent velocities and small-scale virialised mo- 
tions. This will tend to skew the correlatio n functi o n and 
boost the value of ro. The linear analysis of iKaiseJ (Il987l) 
implies a boost in the redshift space correlation function by 
a factor of 1.2 JPeacocklll999ah . Because of uncertainties in 
this effect, and because there are other large uncertainties, 
we have not corrected for this effect. 

Finally, to allow for redshift evolution in the clustering 
of radio galaxies, an evolution term can be incorporated into 
the definition of the correlation length: 



roiz) 



-(3 + £-7)/7 



:ro(0)(l + 2) 

aIl l200.?l : lBro 



(5) 



physical or co-moving co-ordinates respectively. e=7-l rep- 
rese nts growth of clustering under linear perturbation the- 
ory l)Peebles|ll98ff l. Unless otherwise stated, in all following 
analysis, ro(0) refers to the correlation function at 2=0 as- 
suming growth of clustering under linear theory. 



5.3 Calculating the errors on the two-point 
correlation function 

The Poisson errors on the correlation function at a given 
separation r are given by: 



cre(r) 



1 + eM 



(6) 



(e.g. iPeacock fc Nicholson|[l99lh . Poisson statistics assume 
that all pairs are independent of each other which is not the 
case. The Poisson errors will therefore be an underestimate 
of the true error. We also used the bootstrap method to 
estimate the err ors (e.g. Barrow ct al. 1984; Line ct al. 198^; 
iMo et al.llT993) . This involved re-sampling of the data to 
generate further data sets with a population distribution 
identical to that of the real data set. We did this by selecting 
a random sample of A'^ objects from the real catalogue and 
re-calculating the two-point correlation function for the new 
sample. The bootstrap estimate of the standard deviation at 
of the re-sampled correlation function ^* is given by: 



JV-l 

E 

fc=0 



N - 1 



(7) 



fe.g. lOverzier et alJl200.?l : iBrown et alJl200.?f l. e=0 and 6=7- 
3 represent scenarios in which the clustering is fixed in 



where < ^* > is the mean value of ^* in separation bin i. 
We calculated this for 1000 re-sampled datasets (A=1000). 
Because some galaxies were counted more than once, there 
were a large number of real-real pairs in the smallest sepa- 
ration bin. Although this is not therefore a good estimate 
of the mean, it does provide a good estimate of the internal 
error. 

Because the bootstrap method calculates the internal 
variance of the sample, the errors will always be an underes- 
timate of the true error as they don't include the error in the 
underlying number density. However, we have already cor- 
rected for this from the model redshift distribution. The only 
uncertainty will be from the error in our knowledge of the 
underlying number density. We have already estimated the 
la errors associated with the total number of radio galaxies 
expected in each sample (Table. EJ. This was determined 
from the itlcr errors on the model redshift distribution cal- 
culated in Sec. 14.11 To estimate our total error, we added 
the fractional errors in quadrature (including twice the frac- 
tional error on < n > because ^(r) oc< n >^). 



5.4 Estimating ro 

Because each galaxy can contribute to multiple close pairs, 
the correlation function obtained in neighbouring bins may 
be correlated and hence in order to obtain the correlation 
length it is not strictly correct to simply calculate the best 
fit to a power-law function (Equation. |2l . Ideally we would 
obtain the full covariance error matrices using our bootstrap 
error estimate: 
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E 



i=0 



<£.t>m-<c,>) 

N -1 



(8) 



However, to calculate this, ea c h of t he data points must be 
independent. iMadgwick et al.l (|20o3) overcome this problem 
by dividing their survey area into contiguous regions and 
selecting areas at random (instead of the standard boot- 
strap method in which galaxies are selected at random). 
A s imilar approach is that of jackknife re-sampling (see 
e.g. IZehavi et aT]|2002h . Because our sample volume was 
too small to divide into enough regions (each required to 
be larger than the separation distances under considera- 
tion), we instead chose to estimate the value of the correla- 
tion length by using the Levenberg-Marquardt technique to 
perform a non-linear least squares fit of a power-law func- 
tion (Equ ation. |2J to the binned correlation function data 
jPress et al. 1992*1. In this way, we incorporated the cal- 
culated errors rather than assuming simple Poisson errors. 
Because we were using only the diagonal terms of the covari- 
ance matrix, we cannot guarantee that we have obtained an 
unbiased estimate of the correlation length. However, the 
goodness of fit suggests that any such bias is negligible. 



5.5 Relating the correlation length to the bias 

The simplest models relating the galaxy distribution to that 
of the underlying dark matter is a linear bias: pgai — fexpniass 
where 6 is a constant bias factor. The bias is therefore related 
to the correlation function by: 



CgaK'') = b ^mass(r-). 



(9) 



This can be understood intuitively: clusters of galaxies will 
form preferentially at the sites of high-density peaks in the 
primordial density field, and will consequently have a larger 
bias. 

The bias of radio galaxies 6rg with respect to the bias 
of another type of galaxy, fegai is therefore related to their 
scale-lengths ro and power law indices 7 by: 



(ro 



"gal 



(ro gal 



lTgal/2 



(10) 



6 RESULTS 

In this section, we discuss the correlation function and its 
corresponding correlation length obtained for each of the 
TONS sub-samples. We summarise the correlation lengths 
obtained along with values obtained in other studies in Ta- 
ble. □ 



6.1 The TONS08 and TONS12 sub-samples 

Fig. El shows the two-point correlation function calculated 
for the TONS08 and TONS 12 sub-samples. In the case of 
TONS08, there are less radio galaxies than predicted. The 
number of radio galaxies in the random sample is effectively 
increased and the correlation function has therefore been 
corrected downwards. Conversely, the TONS 12 sample has 
more radio galaxies than predicted and the correlation func- 
tion was corrected upwards. 

In both cases, a power-law function provided a poor fit 



to the correlation function. This is due to the presence of 
an excess number of pairs of galaxies with separations of 
30-80 Mpc and is probably caused by the presence of super- 
structures in the sub-samples (see Sec. 16.211 . The error bars 
become large in the low separation bins due to small num- 
ber statistics. The bootstrap errors also become larger when 
the discrepancy between n and < n > becomes large. This 
is due to the multiplication factor of (n/ < n >)^ in equa- 
tion. |1] for each bootstrap realisation. Because of the large 
uncertainties, we fixed 7 to the canonical value of 7=1.8. We 
fitted a power-law function to the correlation function out 
to separations of r-=200 Mpc. Our best fit values for the cor- 
relation length are r-o(z)=8.0±2.4 Mpc and ro(z) = 11.3±1.8 
Mpc (Icr) for the TONS08 and T0NS12 sub-samples re- 
spectively. The best fit power-law models are overplotted in 
Fig. |S| Assuming that the sub-samples have a median red- 
shift of z—0.3 and that clustering grows according to linear 
theory, we used Equation. |3 to determine that the equiva- 
lent correlation length at z—Q are r-o(0) = 10.2±3.0 Mpc and 
7'o(0) = 14.3±2.3 Mpc (Icr). The higher value for the correla- 
tion length in the TONS 12 sample is probably in part due 
to the presence of an excess number of pairs of galaxies with 
separations of 10-20 Mpc which will tend to pull the fitted 
function up to higher values of tq. 

To obtain better number statistics, we computed the 
correlation function for the combined TONS08 and T0NS12 
sub-samples. Fig. |7| shows that the combined correlation 
function has much smaller errors than is obtained for the 
separate sub-samples. This is partly because of the increased 
sample size and also because when TONS08 and TONS 12 
are combined, the number of radio galaxies observed is very 
similar to the number expected (see Table. I^J. 

The best fit correlation length for the combined sam- 
ples is ro(2)=8.7±1.6 Mpc (Icr) (fitting the power-law func- 
tion out to separations of 30 Mpc). This corresponds to 
ro(0) = 11.0±2.0 Mpc (I cr). Using Equ ation. HTI and taking 
ro gai(0)=8.1±0.9 Mpc jNorberg et aU.2002.'l and fogai^l.O, 
this corresponds to a bias of ferg=1.3±0.3. 



6.2 Super-structure regions 

By re-dividing our combined sample into those within and 
out of super-structure regions, we directly determined if 
there are any differences in the clustering properties of ra- 
dio galaxies in different large-scale environments. We com- 
bined the TONS08 and TONS 12 sub-samples and redivided 
the combined sample into radio galaxies that were within 
super-structure regions and those that were not. We de- 
fined super-structure members as radio galaxies with red- 
shifts in the intervals 2=0.233-0.285 and z=0.33-0.367 for 
the TONS08 sub-sample and 0=0.21-0.25 and 2=0.285-0.34 
for the T0NS12 sub-sample. Using this criteria, this re- 
sulted in 50/84 and 42/107 radio galaxies defined as be- 
ing in a super-structure and 34/84 and 65/107 radio galax- 
ies defined as being outside of super-structures in TONS08 
and TONS 12 respectively. The total sample size was 92 ra- 
dio galaxies within super-structures and 99 out of super- 
structures. From our model redshift distribution, we cal- 
culated that we would expect 22.3±2.4 and 24.6±2.4 ra- 
dio galaxies within the super-structure regions and 72.8±3.8 
and 73.0±3.8 radio galaxies outside super-structure regions 
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Survey 


mean z 


ro(z) (Mpc h^i) 


ro{z) (Mpc) 


ro(0) (Mpc) 


b/bopt 


^rbcre ct aL (2002) (L* ellipticals 2dfGRS) 


local 


5.7±0.6 


8.1±0.9 


8.1±0.9 


1.0 


Norbere ct al. ('2002') fbriehtest ellioticals 2dfGRS~) 


~0.1 


9.7±1.2 


13.9±1.7 


15.1±1.9 


1.7±0.3 


^alton ct aL (1994) (APM AROl clusters) 


local 


14.3±2.4 


20.4±3.4 


20.4±3.4 


2.3±0.5 


"croft ct al' (1997') CAPM AR02 clusters') 


local 




3o.4±li:i 


30.4±}i;i 


o o_|_1.9 


Peacock & Nicholson (1991) (local radio ealaxies') 


local 


11. Oil. 2 


15.7dzl.7 


15.7±1.7 


1.8±0.3 


Blake & "Wall (2002") (NVSS radio ealaxies') 


~1.0 






8.6 


1.05 


Maeliocchetti et al. (2004') (2dFRS radio ealaxies) 


~0.1 




13.0±0.9 


14.2±1.0 


1.7±0.3 


TONS08 


~0.3 




8.0±2.4 


10.2±3.0 


1.2±0.4 


TONS12 


~0.3 




11.3±1.8 


14.3±2.3 


1.7±0.4 


TONS08/12 


~0.3 




8.7±1.6 


11.0±2.0 


1.3±0.3 


Lacv (2000) 


~0.3 


17.9±7.4 


12.5±5.2 


15.9±6.6 


1.5±0.7 


TONS08w/16w/Lacv (2000) 


~0.3 




8.8±4.4 


11.2±5.6 


1.4±0.8 



Table 7. Tabic summarising the measured two-point correlation length, tq{z) (h^^ Mpc), the measured two-point correlation length 
assuming h=0.7, ro{z) (Mpc), the correlation function corrected to 2=0 assuming growth of clustering by linear theory, ro(0 ) (Mpc) 
and the corresponding bias with respect to elliptical galaxies in the 2dfGRS iNorberg et al.ll2002^ . iNorberg et al.l l2002^ look at 
clustering properties of different types of galaxies in the 2dfGRS. We take the clustering lengths of L* ellipticals. We assume this is a 
local measurement and use this value as bopt. The brightes t ellipticals are defined as Mb.-51ogio/i-l=-21.0-22.0. To convert to ro(0) 
we assume the median redshift of this population is 2=0.1. iBlake W all 12002) calculate a rough value for the angular clustering of 
NVSS (si. 4 >10 mjy) radio galaxies. The value ro(0) has already been corrected for growth of clustering under linear theory assuming 
a median redshift of 2=1. 




Figure 6. The two-point correlation function for the TONS08 and TONS12 sub-samples. The data are binned in 10 Mpc separation 
bins. The errors plotted are obtained from a combination of the bootstrap method and the error due to the uncertainties in the predicted 
number of radio galaxies. The diamonds show the extent of the Poisson errors. The solid line shows the fitted power-law correlation 
function with ro =8.0 Mpc (TONS08) and ro =11.3 Mpc (TONS12) and a fixed slope of 7=1.8. The dotted Une shows the correlation 
function fitted bv lPeacock fc NicholsonI il99ll) for local radio galaxies evolved to 2=0.3 assuming linear theory (ro=6.1 Mpc). 



in the TONS08 and TONS12 sub-samples respectively (see 
Table. OJ. 

The two-point correlation functions for the combined 
TONS08 and TONS12 super-structure and non-super- 
structure regions are shown in Fig. |H| In these cases, we 
calculated the two-point correlation functions without cor- 
recting for the underlying number density of radio galax- 
ies because we were determining differences in the shape 
and strength of the correlation function in over- and under- 
dense regions rather than trying to determine a universal 
correlation length. The two-point correlation functions have 
very different shapes. The two-point correlation function in 
non-super-structure regions appears to conform to a well- 
behaved power-law relation whereas the correlation function 



in super-structure regions strongly deviates from a simple 
power-law. The enhanced number of pairs at large separa- 
tions seen in the correlation functions of both the TONS08 
and TONS12 samples are due to radio galaxies in super- 
structure regions. This is unsurprising given that we would 
expect a clustering signal on larger scales in regions that we 
see super-structures. WiLd__et al. (2004) use the relative bias 
of galaxies in the 2dFGRS to rule out an exact linear bias. 
The deviation from a simple power-law in super-structure 
regions supports the view that bias may be a more compli- 
cated function. 

The two plots illustrate nicely how a power-law is only 
a good fit either out of super-structure regions or in surveys 
that are less biased and large enough for their contribution 
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Figure 7. The two-point correlation function for the combined TONS08 and TONS12 sub-samples. The data are binned in 10 Mpc 
separation bins. The errors plotted are obtained from a combination of the bootstrap method and the error due to the uncertainties in the 
predicted number of radio galaxies. The diamonds show the extent of the Poisson errors. The solid line shows the f itted power-law corre- 
lation function with ro=8.7 Mpc and a fixed slope of 7=1.8. The dotted line shows the correlation function fitted bv lPeacock fc NicholsonI 
il99 J) for local radio galaxies evolved to 2=0.3 assuming linear theory (ro=6.1 Mpc). 



to be negligible. In a survey consisting of biased tracers of 
the mass, such super-structures will be common and thus a 
power-law will never be a good approximation. 

6.3 The TONS08W, TONS16w and Lacy 
sub-samples 

We also computed the two-point correlation function for the 
sub -samp l es wit h a higher radio flux-density limit. 

iLacvl feOOd) analysed an independent sample of similar 
(si.4 > 20 mjy) radio galaxies in a different area of the sky 
but with a much smaller sample size. To compare our use 
of a different method and in particular, a different redshift 
distribution, we compared the correlati on length that we ob- 
tain with that obtained bv lLacvl i2000l) . Without correcting 
for the fact that the region is under-dense, we obtained a 
correlation length of r pfz) = 15.3 =b 7.1 Mpc. The correla- 
tion length obtained bv lLacvl l|200(tl is ro = ITI^j Mpc. Our 
results agree with this well within the error bars. When we 
corrected the correlation function to that for the underlying 
number density, we found tq^z) — 12.5±5.2 Mpc. This shows 
that the effects of calculating ^ over a small survey volume 
are large. Assuming a median redshift of 2=0.3, this corre- 
sponds to a correlation length at z=0 of ro(0) = 15.9 ± 6.6 
Mpc. 

As demonstrated for the lLacvl (|200nh sample, the small 
numbers combined with the smaller number density of the 
higher flux density radio galaxies in the TONSOSw and 
T0NS16w sub-samples means that the errors are too large 
to determine a reliable value of ro for each individual sub- 
sample. We therefore combined the TONSOSw, TONS16w 



and lLacvl (l2000l) s amples to imp rove the small number statis- 
tics. Although the iLacvl teOOd) sample has slightly different 
selection criterion, the populations are sufficiently similar to 
add the samples together to increase the number statistics 
and still obtain meaningful results. 

The combined two-point correlation function is shown 
in Fig. |n| Because of the possibility of incompletenesses in 
the TONSOSw sub-sample, we assumed that the underlying 
number density of radio galaxies in the TONSOSw survey 
was equal to the number observed (Equation. |1J. The error 
bars are too large to obtain a reasonable estimate of ro. A 
number of reeisons may be accountable for the lack of sig- 
nal and the large error bars. An intrinsic reason may be 
the smaller number density of the higher flux density radio 
galaxies. Perhaps there are too few luminous radio galax- 
ies triggered at small separations. Indeed, there are only 
four, one and zero radio galaxy pairs with separations < 
10 Mpc in the lLacvl tOO(t). TO NS16w and TONSOSw sub- 
samples respectivelv. The lLacvl ll200Q'l sample gives a par- 
ticularly high clustering signal (even after it is corrected 
for being an under-dense region). This is in contrast to the 
TONSOSw and T0NS16w sub-samples which appear to have 
much weaker clustering signals. 

We obtained a clustering length of ro(2)=8.S±4.4 Mpc 
(l(j). Assuming a median redshift of 2=0.3, this corresponds 
to a correlation length at 2:=0 of ro(0) = 11.2±5.6 Mp c (la). 
This is consistent with the flndings of Magliocchct ti et alJ 
(2004) who found no difference in the clustering properties 
of r adio- AGNs of d i fferen t radio luminosity in the 2dFGRS 
and lOverzier et al] i2003l) who found that the angular cor- 
relation function of NVSS radio galaxies is approximately 
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Figure 8. The two-point correlation function for the radio galaxies within (left) and outside of (right) super-structure regions in the 
combined TONS08 and TONS12 sub-samples. The data are binned in 10 Mpc separation bins. The errors plotted are obtained from a 
combination of the bootstrap method and the error due to the uncertainties in the predicte d number of rad io gala xies. The diamonds 
show the extent of the Poisson errors. The dotted line shows the correlation function fitted bv lPeacock fc Nicholson for local radio 

galaxies evolved to z=0.3 assuming linear theory (ro=6.1 Mpc). 



constant between fiux de nsity limits of si.4=3 - 40 mjy. 
iMagliocchetti et al.l (|2004) concluded that once radio activ- 
ity is triggered in low luminosity radio galaxies, there is no 
evidence for a connection between black-hole mass and the 
level of radio output. However our errors are so large that 
one can rule out only very dramatic differences in the clus- 
tering properties of the two samples. 



7 DISCUSSION 

By combining the TONS08 and TONS 12 sub-samples, we 
found a correlation length of ro{z)=8.7±1.6 Mpc. Even af- 
ter correcting for the increase in clustering under linear the- 
ory between z ~0.3 and the local Universe (ro(0) = 11.0±2.0 
Mpc), this value is significantly small er than the local value 
of ro (0)=15.7±1.7 Mpc derive d bv llPeacock fc NicholsonI 
I199J) and that obtained by IMagliocchetti et alJ l|2004fl 
(''o(0) = 14.2±1.0 Mpc). The reason for this discrepancy is 
unclear as we should be looking at a comparable sample in 
terms of the mean radio luminosity of the population (see 
Fig. Iim . This may be evidence for a redshift evolution in 
the clustering of low luminosity radio galaxies but we note 
that the sign of the supposed effect (radio galaxies are less 
biased at z ~0.3 than at z ~0) runs counter to the inferences 
drawn from systematic chan ges in the environ ment of weak 
radio galaxies with redshift llHill fc Lillvlll99ll) . However, a 
number of different effects could contribute to this apparent 
discrepancy including a marginal statistical difference (see 
errors on biases in Table. |7|l, small number statistics, redshift 
space distortion effects and the assumptions that have to be 
adopted when transforming the clustering measurements at 
different redshifts to that in the local Universe. 

We obtained a larger correlation length than that 
obtained from the projection of t he angular correla- 
tion function of NV SS radio galaxies (iBlake fc Walll 120021: 
lOverzier et aljl2003l) . Any discrepancy in this case is easier 
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Figure 9. The two-point correlation fu nction for the combined 
TONSOSw, TONS16W and iLacvl i200(i) samples assuming that 
the predicted number of radio galaxies in TONSOSw is equal to 
the number observed. The data are binned in 10 Mpc separa- 
tion bins. The errors plotted are obtained from a combination 
of the bootstrap method and the error due to the uncertainties 
in the predicted number of radio galaxies. The diamonds show 
the extent of the Poisson errors. The solid line shows the power- 
law correlation function with ro=8.8 Mpc and a fixed slope of 
7=1.8. The dotted line show s the correlation function fitted by 
IPeacock fc Nicholso 3 <199ll) for local radio galaxies evolved to 
2=0.3 assuming linear theory (ro=6.1 Mpc). 



to explain: Fig. ITOl illustrates how the lBlake fc WalJ (|2002h 

sample is dominated by more powerful FRH r adio galaxies 
near the break in the RLF. lMcLure et alJ ll2004ll showed that 
although the scale size and optical luminosity of the host 
galaxy is correlated with the radio power for FRH's, this 
relation does not hold for FRI radio galaxies. If FRI radio 
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galaxies have more luminous host galaxies, they should also 
have a larger clustering length. The radio activity in these 
populations may be triggered by a different mec hanism. This 
is hinted at by the recent result of IClewlev fc Jarvis (200^ 
who found that although it is well known that the co-moving 
number density of powerful FRII radio galaxies falls dra- 
matically at low redshifts, the co-moving number density 
of lower luminosity AGN remains constant. The clustering 
strength of FRII radio galaxies may therefore be more com- 
plex than simply tracing the clustering of the underlying 
host galaxy population. However, these measurements may 
have significantly underestimated the correlation length for 
two reasons. Firstly, deep cleaning of the NVSS images may 
have caused a deficit of galaxy pairs with separations 6 --^0.2 
deg, which may have resulted in a spuriously low p ower- 
law fit and hence correlation length (iBlake et alJl20o3l . Sec- 
ondly, the assumed redshift distribution used to deproject 
the angular correlation function has an unphysical spike at 
low redshifts which would also have suppressed the correla- 
tion signal. When these effects are accounted for, the correla- 
tion length obtained was ro(0) = 15.1±0.7 (C. Blake, private 
com munication) . 

iNorberg et al.l ll200ll) showed that the strength of 
galaxy clustering depends on the galaxies intrinsic optical 
luminosity. The mean luminosity of the TONS08/12 ra- 
dio galaxi es is ^^2.5 L<,. From the relation: b/b* = 0.85 -I- 
0.15L/L* iNorberg et all i200ll) . we expect the sample to 
have a relative bias of b=1.2. From INorberg et al.l ||2002|) 
fig. 10, the expected bias of 2.5L* early-type galaxies has 
a slightly higher bias of 6=1.3. Using equation. 10, the 
TONS08/12 (linear theory evolution corrected) clustering 
length corresponds to a bias of &=1.3±0.3. Although the 
uncertainties are large, this suggests that the clustering 
strength of the TONS radio galaxies is consistent with that 
of the underlyin g host galaxy popula tion. If, as suggested 
by the results of McLurc ct al. '2004j , the clustering signa l 
measured bv >Blake fc Wall 12002.) and lOverzier et ali(l2003l) 
was dominated by galaxies with smaller (~ L*) optical lumi- 
nosities, we would indeed expect them to measure a smaller 
clus tering signal. 

iBrand et alj (2003) postulated that the radio galaxy 
bias of low luminosity radio galaxies needs to be h ~3-4 
to explain the number of super-structures found within the 
TONS08 volume. In this paper, we have shown that redshift 
spikes are a common phenomena in radio galaxies redshift 
surveys at NVSS flux density levels and z ~0.3 (correspond- 
ing to FRI-like radio galaxies well below the RLF break). 
The d iscrepancy in the bias inferred here and bv lBrand et alJ 
i2003h may be related to the inadequacy of a simple linear 
bias model for the way radio galaxies trace the underlying 
dark matter. 



8 CONCLUSIONS 

The results of this investigation can be summarized as fol- 
lows. 

• We have demonstrated that redshift spikes correspond- 
ing to super-structure size overdensities are a common phe- 
nomenon within radio galaxy redshift surveys. 

• We computed a universal redshift distribution for radio 



galaxy redshift surveys of given optical magnitude and radio 
fiux density limits. 

• The spatial clustering of low luminosity (predominantly 
FRI) radio galaxies at z ~0-0.5 has been measured directly 
for the first time. 

• We used our model redshift distribution to correct the 
clustering signal for our sample being under- or over-dense 
due to large-scale structure. 

• The best fit correlation length for our low fiux density 
limit (si.4 >3 mjy) sample is ro(2)=8.7±1.6 Mpc, corre- 
sponding to ro(0) = 11.0±2.0 Mpc at z=Q (assuming an in- 
crease in clustering under linear theory). This is consistent 
with the clustering strength of the underlying host galaxy 
population. 

• We found that the clustering of low luminosity radio 
galaxies is not significantly different for our low and high 
fiux density limited samples. 

• By comparing the two-point correlation function for ra- 
dio galaxies within and outside of super-structure regions, 
we showed that a power-law function is only a good approxi- 
mation outside super-structure regions or in surveys that are 
less biased and large enough for their contribution to be neg- 
ligible. In a survey consisting of biased tracers of the mass, 
such super-structures will be common and thus a power-law 
will never be a good approximation. 
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